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Abstract 

Using closed-time path two-particle irreducible coarse-grained effective action (CTP 2PI CGEA) 
techniques, we study the response of an open interacting electronic system to time-dependent 
external electromagnetic fields. We show that the CTP 2PI CGEA is invariant under a simultaneous 
gauge transformation of the external field and the full Schwinger-Keldysh propagator, and that this 
property holds even when the loop expansion of the CTP 2PI CGEA is truncated at arbitrary order. 
The effective action approach provides a systematic way of calculating the propagator and response 
functions of the system, via the Schwinger-Dyson equation and the Bethe-Salpeter equations, 
respectively. We show that, due to the invariance of the CTP 2PI CGEA under external gauge 
transformations, the response functions calculated from it satisfy the Ward-Takahashi hierarchy, 
thus warranting the conservation of the electronic current beyond the expectation value level. We 
also clarify the connection between nonlinear response theory and the WT hierarchy, and discuss 
an example of an ad hoc approximation that violate it. These findings may be useful in the study 
of current fiuctuations in correlated electronic pumping devices. 
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I. INTRODUCTION 



In this paper, we are interested in the response of an open electronic system to time- 
dependent external electromagnetic fields. Although our analysis is quite general and can 
be applied to a variety of condensed matter systems, we have in mind a system of strongly 
interacting electrons subjected to external driving fields and in contact with two ideal reser- 
voirs of noninteracting electrons. This could be a picture of the so-called interacting elec- 
tron pumping devices, which are currently attracting much interest both experimentally 
and theoretically^. In these devices, a direct current can be generated by applying slowly- 
oscillating external fields to the central electrons, even in the absence of a bias voltage 
difference between reservoirs. These are very interesting systems to study from the theo- 
retical point of view, since they represent a unique and challenging combination of strongly 
correlated particles and quantum transport, and for which the study of current conservation 
beyond the expectation value level is nontrivial. A typical situation where one needs to go 
beyond the expectation value of the current is when dealing with current fluctuations. It is 
with these applications in mind that the following considerations were developed. 

The main purpose of this work is to determine the basic requirements that a fleld- 
theoretical approach to open driven systems must satisfy in order to produce current conserv- 
ing results (in the sense of the WT hierarchy developed in Sec. |Tl]), in transport calculations 
going beyond linear response. Another aim we have in mind is to clarify the close relation 
that exists between current conservation and response theory, especially in the nonlinear 
regime. In order to analyze these issues, we combine the so-called external gauge invariance 
method with the closed-time path 2PI coarse-grained effective action (CTP 2PI CGEA), 
suitable for the description of strongly interacting quantum open systems both in and out 
of equilibrium^>^>^>^i^. 

The area of nonequilibrium physics in interacting systems is gaining increasing interest 
nowadays, particularly time- dependent quantum transport in correlated system3^'^'^»^»^. 
The so-called time propagation method! constitues a signiflcant advance in this area. It con- 
sists in flrst determining the (interacting) equilibrium Green function and then propagating 
it by using the Kadanoff-Baym equations. This is equivalent^i^ to solving the Bethe-Salpeter 
equation for the particle-hole propagator (which is related to the flrst order response function 
of the system), but numerically much less expensive. This is a powerful method in which 
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external fields are treated exactly to all orders while many-body interactions are treated 
perturbatively. In this paper, we adopt another method to study nonequilibrium transport 
through correlated systems, in which both the external fields and the many-body interac- 
tions are treated perturbatively. Since it is based on an expansion of the Green function in 
powers of the external fields, it is only valid for weak external fields. On the other hand, it 
allows us to focus on current conservation beyond the expectation value level (i.e. conserva- 
tion of current fiuctuations). In this aspect the present work goes beyond previous analysis 
such as those described in Refs. 

The paper is organized as follows. In Sec. [Ill we present the Ward-Takahashi (WT) 
hierarchy, which is the most general form of current conservation beyond the expectation 
value level, for a system driven by external fields. We discuss two possible ways (based on 
effective action techniques) of generating approximations to the non-equilibrium many-body 
problem that satisfy the WT hierarchy. One of such methods, relying on the external gauge 
invariance of the effective action, is used throughout the paper. In Sec. Hn] we introduce 
the 2PI CGEA of the system, closely following Ref. 3. In Sec. IIVI we prove that the exact 
and truncated effective actions are external gauge invariant. Using these theoretical tools, 
we show in Sec. |V]that for any approximation to the Schwinger- Dyson equation (obtained 
from a truncation of the loop expansion of the 2PI CGEA) there exists a corresponding 
approximation to the Bethe-Salpeter equations (which give the vertex functions), such that 
the WT hierarchy holds. The WT identities are systematically obtained from the external 
gauge symmetry of the 2PI CGEA. In this section we also clarify the relation between 
the WT hierarchy and current-conserving nonlinear response theory. In Sec. |Vl] we show, 
using a simple example, how ad hoc approximations to the 2PI CGEA (not resulting from 
a truncation of its loop expansion) violate the WT hierarchy. A brief summary is given in 
Sec. IVm 



II. WARD-TAKAHASHI HIERARCHY AND CURRENT CONSERVATION 

In many-particle systems and particularly in quantum transport theories, n-point vertex 
functions play a fundamental role, since they represent generalized currents which satisfy 
the hierarchy of Ward-Takahashi identities'^. This hierarchy is satisfied to all orders in the 
exact theory, thus guaranteeing local gauge invariance and the conservation of the associated 
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charges. 

There are two equivalent ways of obtaining the WT identities. The first one, due to 
Rivier and Pelkap^ (See also Ref. Iisl ). relies on the equation of motion of the n-particle 
Green's function (GF). In the second one, due to Kadanoff and Baym^"^'^^, the two-point 
GF is expanded in powers of a fictitious non-local external field (for a recent development 



see Ref. 



101 ). To zeroth order in the expansion, the classical continuity equation is recovered 



(i.e. the mean current is divergenceless) ; to first order in the external field, the usual WT 
identity relating the three-point vertex to the two-point propagator is obtained. Higher 
order terms in the expansion reproduce the hierarchy of WT identities obtained by the first 
method. We note that the second approach is not directly related to Baym's $— derivable 
approximation^ii^'i^, because in principle no approximation is involved. It is just an equiva- 
lent way of deriving the WT hierarchy which is based on a functional approach and not on 
the equations of motion for the propagators. 

In its most general form, the generalized continuity equation (WT hierarchy) can be 
written as (we use the simplified notation 1 = (ti,ri) and employ Schwinger-Keldysh non- 
equilibrium iormalismr^'^'^'^'^'^) : 

d;Al^il...n,l'...n';z) = 

V^ei-d^z - n')Gn{l ... n, 1' ... (n - l')z) + ... 

+ {-ir5{z-l')Gn{l...n,2' ...z) + ... (1) 
+ 5{z - n)G„(l ...{n- l)z, I' ...n')- ... 
-{-ir5{z-l)Gn{2...z, l'...n')} 



(2) 



where 

Af„)(12...n,l'2'...n';^) = 

< T,/(z)^(l)^(2) . . . ,p{n)^\n') . . . 'iP\2')ij\l') > 
is the (n -|- l)-point vertex function with current insertion at z = (tz, r^), 
f{z) = -e lim D^'iz, z')^\z')^{z) ; 

D\z,z') = i2t)-\Vi-Vi,) /i = 2 = 1,2,3 , (3) 

D°{z,z') = l fx = 
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and Gn are real-time propagators defined as usual 



...n,l'...n) = 

r(") < r,V^(l) . . . ^{n)ilj\n) . . . > 



(4) 



Note that A^f^^^Jz) =< jt'[z)>, and that 



D\z,z')[A] 



i{Vl-Vl.-ie[A\z) + A\z')]}/2 



(5) 



in the presence of external fields. 

In Eqs. ^ and (jlD, is the closed-time path time ordering operator, ip and ip"^ are field 
operators in the Heisenberg representation, and < . . . > stands for an average taken with 
respect to the density matrix of the system in the remote past. The classical continuity 
equation and the usual WT identity^i correspond to the cases n=0 and n=l of Eq. ([1]), 
respectively: 



being G'(l, 1') the single-particle propagator. 

As already mentioned, in the exact theory Eq. ([T]) holds for arbitrary n. The hierar- 
chy of n-particle propagators take into account that field excitations can be created and 
annihilated in the second-quantization formalism, thus acting as source terms for the gen- 
eralized currents represented by the vertex functions. In other words, in the exact theory, 
particle number is strongly conserved, not only in the mean. This fact has an important 
consequence in quantum transport theories based on response theory to an external field, 
which is the main subject of this work. For a many-particle system in the presence of an 
external field, the hierarchy for the exact propagators implies that the current, as defined by 
Eq. ([3]), is conserved to all orders in the external perturbation, with linear response theory 
corresponding to the special case n=l. 

It is impossible, in general, to obtain the propagators of an interacting field theory exactly, 
and some approximations must be made. In the weak-coupling regime, no conflict arises be- 
tween conservation laws and perturbative expansions based on Feynman diagrams, because 
conserved quantities are conserved order by order in the expansion^^. On the other hand, 
for a strong coupled theory re-summation concepts are usually needed^i^^, and warranting 




(6) 



and 




(7) 
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conservation laws then becomes a nontrivial issue. Symmetries satisfied at the exact level 
may not be satisfied by approximate propagators, thereby violating (generalized) current 
conservation dictated by the hierarchy in Eq. ([T]). 

A systematic way of generating conserving approximations (at the classical level, i.e. 
n=0) was given by Baym^ and corresponds to his well-known ^-derivable scheme. The 
self-energy is obtained from a functional $ consisting of an infinite series of two-particle 
irreducible (2PI) closed diagrams, constructed from full propagators and bare vertices^. 
The solutions obtained from truncating the $ functional are such that, if $ is invariant 
under a simultaneous symmetry transformation of the classical field and propagator, the 
expectation values of the respective Noether currents are conserved^i^i2^. We emphasize that 
this situation corresponds to the case n=0 of the WT hierarchy, given by Eq. Therefore, 
the conservation of generalized currents encoded in the WT hierarchy is not automatically 
warranted in this approach. Moreover, in certain situations (e.g. in theories exhibiting 
spontaneous symmetry breaking) the hierarchy is already violated at the level of the self- 
energy (corresponding to n =1) which does not satisfy the Nambu-Goldstone theorent^'^'^. 
The reason is simple: the n-point functions obtained by functional differentiation in a 
derivable approximation may not be equal to the one-particle irreducible^ (IPI) functions 
that satisfy the WT hierarchy. 

To our knowledge, there are two possible approaches to overcome this problem present 
in ^-derivable approaches^. In the one put forward by van Hees and KnolP^, a non- 
perturbative approximation for the IPI quantum effective action is obtained on top of a 
self-consistent solution to Schwinger- Dyson equations (SDE) derived from the truncated 
2PI effective action. The vertex functions obtained in the usual way from this IPI effective 
action fulfill the WT identities. The extra terms not accounted for in the ^-derivable 
approximation and required to satisfy the identities are encoded in a Bethe-Salpeter equation 
(BSE) and higher vertex equations (see Ref. l27| for a calculation of conductivity in QED). 



The second approacb^"^ is based on the concept of external gauge invariance, and provides 
a systematic way of generating consistent SDE and BSEs that automatically satisfy WT 
identities. Most importantly in the context of transport theories of strongly interacting 
systems, the derivation of the SDE and the BSEs can be done, in principle, to any order 
in the external field coupled to the system. Therefore, the latter is especially suited for the 
study of response theory beyond first order and its relation to current conservation beyond 
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the expectation value level, and will be employed here. 

III. CTP 2PI COARSE-GRAINED EFFECTIVE ACTION 

In this section, we introduce the closed-time path 2PI effective actionSi'^i^'^'^ii^ for a spe- 
cific fermionic system open to an environment. A detailed presentation of effective action 
techiques and the CTP formalism can be found in Ref. [3. 

The system we are interested in consists of a central region with interacting electrons 
described by fields ip) and coupled to an external field. This central region is connected 
to two ideal reservoirs of noninteracting electrons described by fields (0]j,,0q,), with a = 
[L, R) denoting left or right reservoir. The reservoirs are assumed to remain in equilibrium 
at all times. 

The CTP classical action of the system is S[iIj, tp, 0, cp] = + + Sc, with: 

= CAB^^B-/j,<P'' (8) 

Sc = cab'4^^Tab<P^ + h. c. 

where 

being Uabcd the completely antisymmetrized bare interaction local vertex. For the moment, 
we shall not specify the structure of the four-fermion vertex Uabcd since its precise form 
will not be needed until Sec. I VII Tab is a local coupling parameter between the central 
region and the reservoirs. 

We are using a DeWitt notation^i^ with A = (x,a), x = {tx,rx,(j) and a = (1,2) [or 
(+, — )] being CTP indices indicating the branch within the closed-time contour. For the 
fields describing electrons inside the reservoirs, an additional index a must be included in the 
CTP indices {A, B), but for simplicity we leave it implicit. Repeated indices are assumed to 
be integrated or summed, {ip, ip) and (0, 0) are Grassmann variables and Cab is a CTP metric 
Cab = diag(l,— 1), while cab = CabS{xA,XB)- Cab and Bab are the free CTP propagators 
corresponding to and S^. They satisfy (we set h,m,c = 1 in what follows) 

cabOBbc = i5A 

(10) 

cab^Cbc = 
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where □(!) = idt^ — /?.o(l) while □(!) = idt^ — /io(l), being and /iq single-particle Hamil- 
tonians corresponding to the reservoirs and the central region, respectively. Explicitly, they 
read 

/io(l) = and 

1^ (11) 

ho{l) = -HVi + eA(l)]2 + eAo(l) 

where A^{z) is the external classical field and i = (1,2,3). Note that it is not necessary to 
add CTP indices to the external field, since this is a physical source, and in any case we 
would obtain i = ^4^2- It is important to remark that, since is external, it is not 
treated ClS db dynamical field. 

From the classical action of Eq. ([8]), we form the CTP generating functional Zctp with 
local and bi- local external sources J a and Kab, respectively. Since we are interested in the 
dynamics inside the central region coupled to reservoirs, only the fields ip) describing 
the "system" will be coupled to the external sources. Zqtp contains all the information 
about the non-equilibrium many-body system since all closed-time path correlators can be 
calculated from it^. When written as a path integral over field configurations it reads 

ZcTp[J.K] = J [D^][D<l)][D^][D^]expt{S[^,^,^,<j)] + JAij'' + ij''JA + lKAB^''^''^ 

(12) 

where ( J^, Ja) and Kab are fermionic and bosonic sources, respectively. The measure cor- 
responding to each field, for example ip, actually stands for [-D^/'a] with a = (1, 2) the branch 
index. The CTP boundary condition of the path integral giving Zctp, namely the conti- 
nuity of field histories in the remote future, is implicit. We assume that the initial state is 
prepared in the remote past, corresponding to the m-vacuum2.i2,. Note that {Ja, Ja) act as 
Lagrange multipliers constraining the deviations of ip) from their mean values, while Kab 
constrains their fiuctuations^. This is the main idea behind nPI effective action techniques, 
namely, to treat 1 . . . n— point correlators on the same footing^i^iS^. By using the 2PI CGEA, 
only the one and two-point correlators are treated as dynamical variables and therefore ob- 
tained variationally; higher correlations (with n > 2) are calculated from them. Note that 
the use of the 2PI CGEA to obtain equations of motion for mean fields and propagators 
represents an enormous simplification of the usual perturbation expansion. Moreover, non- 
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perturbative approximation schemes based on expansions of the 2PI EA provide powerful 
calculational tools where standard expansion schemes break down^i^. 

Any correlation function of the many-body system can be obtained from Zctp by func- 
tional differentiation with respect to {J a-, J a) and then setting all sources to zero. For 
example, the CTP two-point propagator of the central region is 

[D4>] [D<P] [I^V] [I^VIV'^^'' exp tP, ^A]) , 

with 6ii^L denoting right and left differentiation, respectively. Note that there is no need 
to use contour-ordering operators in the second line of Eq. f[T^ . since the path integral 
automatically arranges operators in the correct order. 

Since the reservoirs' fields enter the action quadratically, we can integrate them out 
exactly in Zctp, thus defining a new generating functional for CTP propagators belonging 
to the system. This coarse-grained generating functional reads: 

ZcTp[J, = j expi{S[ip, ip] + sources) (14) 

where the effective classical action is given by 

S[^, ^]=S^ = CABi'^Df^i;'' + SrntH^, ^j] (15) 

with 

DaI = Cab + ^^^,ab • (16) 

The self-energy describes the influence of the reservoirs on the particles' dynamics inside 
the system, and it is given explicitly by 



kAD 



Tab^bcTcd (17) 



where B is the given by Eq. (fTOj) . It is a complex quantity whose imaginary part represents 
the tunneling rate of particles from the central region to the reservoirs. Note that can 
be easily calculated since it depends on the equilibrium CTP propagators of the reservoir 
composed of non-interacting electrons. The CTP propagators corresponding to the system 
are calculated from Zqtp using the analog of Eq. (IT^ with the replacement S S. 

From Zctp we define the CTP generating functional of connected propagators Wctp in 
the usual way, i.e. Zctp = expiW, and then apply a double Legendre transform to obtain 



the 2PI-CTP effective action of the system. The sources {J,K) are connected to the mean 
field and full propagators (denoted by G) through 

dW 



SJ 



= 



A 



^ = f^ = (18) 



A 



6W 1 

T^'^AB 



5Kab 2 

where we have used that the fermionic mean-fields (i^yip) = since no symmetry breaking 
occurs. By performing the double Legendre transform on W and then using the background 
field method we can write the CTP 2PI CGEA of the system as^ (Tr and In operations are 
understood in a functional sense) 

r [G, A] = iTi InG - iD-^lG^^ + Pa [G] (19) 

where A is the external classical field introduced before. We note that the CTP 2PI CGEA 
depends explicitly on the external field only through this will be important in what 

follows. 

In Eq. f|T9|) . r2[G'] encodes all quantum corrections and consists of vacuum 2PI closed 
diagrams with full propagators in internal lines and vertices corresponding to a theory with 
shifted classical action + il),ip + (neglecting constant and linear terms), where ip) 
denote fluctuations. The diagrams are vacuum because the mean value of the fluctuation 
fields is zero by construction. Note that because of the vanishing of the mean fields, the 
shifted classical action vanishes when evaluated at {tp, ip). Therefore, the vertices contribut- 
ing to r2 are identical to those of the original classical action, namely, the quartic one with 
U as coupling parameter. The case of noninteracting electrons corresponds in this scheme 
to T2 = 0. 

The Schwinger-Dyson equations for the CTP propagators follow directly from r[G, A] 

-Ikab (20) 



SGab 2 

where the physical case corresponds to vanishing external sources K = 0. Defining the 
self-energy 

Eab = -25r2/5G^^ (21) 
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we can rewrite the SDE in the usual way 



GaI = Dab + ^^AB , (22) 

being E^s one-particle irreducible by construction. 

In this paper, we shall not go beyond second order in the interaction parameter U, 
corresponding in this theory to a three-loop expansion. To this order, r2 reads in our 
compact notation: 

r2 = — -^UabcdGabGcd + 

1^ (23) 

i — UABCoUEFGHGAEGBpGcGGDH 

The first term corresponds to the so-called double-bubble diagram and the second term to 
the basketball diagram. In the SDE, the first term yields the time-dependent Hartree-Fock 
approximation, while the second is non-local and complex (therefore including fluctuation 
damping)^|^>^. We will return to these approximations when proving the external gauge 
invariance of the CTP 2PI CGEA, in Sec. HVl and also in Sec. |Vll 

IV. EXTERNAL GAUGE INVARIANCE OF THE CTP 2PI EA 

The crucial observation that allows us to relate the CTP 2PI CGEA of the open system 
to nonlinear transport through it is that r[G', A] is invariant under a gauge transformation 



of the external field A^. Following Bando, Harada and Kugo (Ref. l28l ). we will call this 
external gauge invariance (EGI) of the 2PI EA. We will first give a brief proof of the EGI of 
the exact 2PI CGEA, where exact means that r2 is expanded to all loop orders in Eq. fll9l) . 
Then, we will prove that even when the loop expansion of r2 is truncated at certain order, 
the truncated 2PI EA still remains EGI. 

Under a local transformation f/(l) = exp ie(f{l), the external field and the full propagator 
transform asr^^ (we omit the CTP indices for the moment) 

A,^A'^ = UA,U-'-t{d,U)U-' 

(24) 

G(l,2)-.G'(l,2) = f/(l)G(l,2)f/-i(2) . 

It is rather straightforward to prove that, if we retain all terms in the loop expansion of r2, 
then the 2PI CGEA is EGI-'^'^^. This can be simply understood recalling that the exact 
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2PI CGEA is precisely the generating functional of IPI propagators. The proof of the EGI 
of the exact 2PI CGEA is based on the fact that, under a gauge transformation of the fields 
{Afj^jip, ip), the CTP generating functional Zctp[J, K] as given in Eq. fll4p is invariant. This 
is because the transformation is equivalent to a change of integration variables in the path 
integral, whose Jacobian is trivial due to the unitarity of the transformation (det U = 1). 
Since the classical (effective) action S[ip,ip] is external gauge invariant, only the source 
terms are transformed. For a local infinitesimal variation of the fields ( = (A^, ■?/'), given 
in compact notation by 

6C = e{0,-tip^,tip^) + {df,ip,0,0) (25) 

we obtain from Eq. f[T^ 



5ZcTP = =< JaS^P^ + 5i)''jA+ 



(26) 



with 

<...>= / [Dij][Dij]{. . .) expi{S[ij,ij] + sources) . (27) 



Using that 



- -j^ + -Kab4^'' 



5i)^ 2' 

J A - Ikab^P"" (28) 



Sip^ 2' 
--Kab 



5 Gab 2 

which follows since W and F are Legendre transforms of each other, we can turn Eq. (126|1 
into an equation for F: 

E ^C.|^ + T^^G- = (29) 
„tt3 ^C" 5Gab 

where the variation of the full propagator is given by 

= S< ij^ij"" >=< Sij'^ij'' + ^^5^^ > . (30) 



According to Eq. ( I25l) . Eq. (13011 is the infinitesimal version of the transformation rule for G 
given in Eq. ( 12^ . Therefore, Eq. (129|) shows that the exact 2PI CGEA is gauge invariant 
under Eq. 
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We turn now to the invariance of the 2PI effective action which results from a truncation 
in the loop expansion of the quantum correction For completeness, we will also show 
explicitly that the first two terms in Eq. fll9l) . the one- loop terms, are EGI. We start with 
the first term of T[G, A], which is clearly EGI: 

Tr InG ^ Tr \n{UGU-^) = Ti[U{\nG)U-^] = Tr InG. (31) 

The EGI of the second term of Eq. ( IT9l) is proved by noting that D^^ transforms as G 
under Eq. (1241) . This follows directly from the equation of motion for the system's free 
propagator C, given in Eq. (fTOj) . and the fact that E,^ is invariant under the external gauge 
transformation. The latter is a consequence of the equation of motion satisfied by B, the 
reservoir's free propagator, and the dependence of S,^ on this propagator, given in Eq. (fTTl) . 

The third term contributing to r[G', A] is r2[G'], the sum of 2PI vacuum diagrams with 
G in internal lines. The proof of external gauge invariance proceeds as before, and relies on 
the transformation law for propagators and the structure of the 2PI EA. To be explicit, we 
rewrite the first term of Eq. fl23l) taking into account that the interaction vertex Uabcd is 
local. Using the locality of the bare vertex we obtain for the double-bubble diagram 

r2^'* = —-^GabGba (32) 

o 

which shows that, as already mentioned, the self-energy derived from it is local. It is clear 
from the above expression that, under the external gauge transformation given in Eq. fl24l) . 
r2^^ is invariant. The structure of higher terms in the loop expansion of r2 is such that the 
EGI holds to arbitrary order. 

The main conclusion of this discussion is that, even though the system is open to reser- 
voirs, its 2PI coarse-grained effective action is, order by order in the loop expansion, invariant 
under the gauge transformation given in Eq. fl2^ . As we shall see, this makes the combi- 
nation of the CTP 2PI-CG effective action approach and the external gauge invariance a 
powerful technique to study strongly interacting driven systems beyond linear response. We 
should remark that in the presence of non-vanishing mean fields, the 2PI CGEA may not 
be EGI order by order in a loop expansion^'^'^. 
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V. WT HIERARCHY FROM EXTERNAL GAUGE INVARIANCE 



We have seen that the CTP 2PI CGEA describing the system is left invariant under 
an external gauge transformation, and that this property holds to all orders in the (loop) 
expansion of the quantum corrections given by This implies that the SDE derived 
functionally from the 2PI effective action is automatically external gauge covariant. Most 
importantly, this is true for an arbitrary external classical field. This will be the key property 
to relate this symmetry with nonlinear response, thus generating approximations which 
satisfy the WT hierarchy to arbitrary order. We will show that the external gauge covariance 
is also inherited by the BSE for vertex functions, which as a consequence of EGI will satisfy 
the Ward-Takahashi hierarchy. 

Gauge invariance of the effective action T[G, A] implies 

r[G,A] = r[G',A'] (33) 

where a prime denotes external gauge transformated quantities according to Eq. (^^. The 
SDE for the CTP propagators is given by ST/SGab = 0, and due to the EGI of the 2PI 
CGEA, it is (external) gauge covariant: 

6T[G,A] _^^_, ST[G',A'] ^^ _ ^^^^ 
6 Gab ^^'ab 

This property of the 2PI CGEA has an important consequence. It implies that the 
variational procedure by which one obtains the SDE from the 2PI CGEA is independent of 
the external gauge. That is to say, if G^^f^] is a solution of the SDE in the background 
classical field A, i.e. 

^ ^0 , (35) 

OGab Gab=Gab[A\ 

then the solution corresponding to a transformed external field A' = UAU~^ — i{dU)U~^ is 
precisely G' = UGU-^: 



Sr[G',A'] 



^G'ab G'^g=UGAB[A]U- 

More explicitly, EGI of the 2PI CGEA implies 



. (36) 



G[A'] = UG[A]U-^ . (37) 

We will now make a connection between EGI and nonlinear response. For notational 

ution G[A] to the SDE can be 



simplicity, CTP indices are omitted in the following. The so 



expanded in powers of the external field (see Refs. 

14 



18 



2G 



3l|): 



G[A]{X, Y) = G'[0](X, Y) + lA^W^ + ^A^A.nf + ^A^A^A^nf^ . . . 

■n !• (38) 
= J2- dl... dnA,,{l) . . . A^„(n)nf-2^;(1, . . . , X, F) , 

n=0 ■ 

where in the second hne we have made exphcit the internal (integration) (1, . . . , ra) and the 
external (X, Y) variables. The dummy indices fi, u, p, etc. of the first line are denoted by 
/ii, . . . , /i„ in the second. 

The "response" functions n(n+2) encode the variation of the full propagator with the ex- 
ternal field. We note that it is possible^^i^, in principle, to expand higher order correlation 
functions similarly to the two-point function as given in Eq. fl38l) . but in this paper we shall 
only deal with G[yl](X, F). Eq. ( |38i) is given in the so-called™-'^ single-time representation 
of CTP correlators, which is simpler for calculations than the "physical" or r/a represen- 
tation. Both representations are related by a similarity transformation characterized by an 
orthogonal matrix— Q = {1 — iO'i)l\pl with 





and are completely equivalent to each other: 

. . .n) =1 ^ Qi-^^cti ■ ■ ■ QinCen-^cti...a„i^ ■ ■ ■ 

^a,...a.(l ...n) = 2'--/'Ql^^ . . . Q^„,„E,.„,J1 . . . 



(40) 



where E{l...n) denotes a n— point CTP function, aj {ij) corresponds to the single-time 
(physical) representation and repeated indices are assumed to be summed. To give an 
example of Eq. (HOj) , for the propagators in the physical 




G= I „ „ I (41) 

and in the single-time 

(5=1"-"- I (42) 




15 



representations we have^: 

Gil = 2^^++ ~^ ^ — ~ ~ '^-+) = 
Gl2 = Ga = -(G'++ — G — y + G^ — — G — ) 

G21 = Gr = - G+_ + G_+ - (43) 

G22 = Gc = + G — + G^ — + G — 1_) 

= G-\ — -|- G — |- = -l- G — 

The first equality in Eq. fH5]) is valid in general: any CTP function in the physical rep- 
resentation with all its indices set to "1" satisfies = 0. The last line results from 

the identity G+- + = + G , which follows from the normalization of the step 

function 9{1, 2) + 9{2, 1) = 1. In these expressions, Gj., Ga and Gc are retarded, advanced 
and correlation functions, respectively: 

G.(l,2) = -^^(l,2)<Ml),^t(2)}> 

a,(l,2) = -ie{2A) < {V'(l),^t(2)} > (44) 
Gc{l,2) = -t<{^{l),ij\2)}> 

where {, } is the anticonmutator and 6 the step function, while 6"++, ,h are the chrono- 
logical (Feynman), antichronological (Dyson), lesser and greater correlations: 

G'++(l,2) = -z<TV'(l)V'^(2)> 
G+_(l,2)=z<^t(2)^(l)> 

(45) 

G'_+(l,2) = -i< ^(1)V^^(2) > 
G__(l,2) = -2<f^(l)^t(2)> , 

with {T,T) denoting the chronological and antichronological time ordering operator. Using 
this transformation rule, it is a simple matter to re-ex pre ss Eq. ( 138|) (and what follows from 



it) in the physical representation (see Refs. 



18 



and 



20I ). but for clarity we will continue 



employing the single-time representation. 

We note that, as discussed in Ref. Q, in the physical representation the observables 
are given by retaining only the "2" component of the n— point function G{1 . . .n), i.e. by 
the correlation functions G2...2- For example, the two-point observable (n = 2) is the fully 
symmetrized correlation function G22{X,Y) = Gc{X,Y), given in Eq. (H3|) . Using the 

16 



relation between the physical and single-time representations, the response of the correlation 
two-point function G22 is immediately obtained 

G22[A]iX,Y) =G22mX,Y) + 

(46) 

zn^2i(x,r;i)A^(i) + ... 

It is worth remarking that the use of the closed-time path method automatically ensures 
that the response functions are causal (see, for instance, Refs. yandlisl}^. 
Returning now to Eq. (1381) . the response functions are given by 

Ul^-;^;il,...,n;X,Y) = 

6G{X, ¥■ A) (47) 
(1) . . . (n) [A^^ ...A^„]=o 

and correspond to [n + 2)— point functions with n current vertices inserted at locations 
(1, . . . , n) where interactions between the current and the external classical field take place. 
In more detail, the response functions for n = 1,2 read 

n^(X, Y; z) = ~i< T,f{z)^{X)^\Y) > and 

(48) 

UT{X,Y;z,w) = -< T,nz)f{w)^{X)^\Y) > 

Note that —iIl2{X,Y) corresponds to the propagator in the absence of the external field, 
which can be written as G'[0](Xj, Yf, Xq — Yq) due to time-translation invariance of the system 
in equilibrium. 

The structure of n(„_|_2) as given in Eq. fHSj) follows because functional differentiation with 
respect to yields a current insertion to which the external field couples. In functional 
language we have that 



6G{X,Y-A) 



6A,,{1)...6A,M 



[D^][d^]^/j{X)^lj{Y) exp2(5[^/',V] + sources) 

Eq. (149!) automatically leads to Eq. (148|1 since 

S 



expi{S[ip,ip] + sources) 



SA,,{1)...6A,M 
= (-0"Jmi(1) • • •iMn('^)exp2(S'[V^,V', A = 0] sources 
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The functions n(„_|_2) are obtained from the SDE by functional differentiation with respect 
to Afj_ (and then setting A = 0). This results in the BSEs for the response functions. For 
example, the BSE that determines lis reads 

^1— = ^1— + ^T1- (51) 
SA^ 6A^ 6A^ ^ ' 

where both sides of the equation are evaluated at A^ = 0. Using that G(1,2)G~^(2, 1') = 
(5(1, 1') (where the coordinate 2 is integrated), which implies 



and defining 



^^-/.4.5G-.,M)^G-.,5,2, (52) 



the BSE for lis can be written as 

^^-^-^ ■ 

Since the self-energy E is defined through it is a functional of the full propagator G and 
therefore the last term in Eq. (1541) involves II^. Namely, 



5A^ 5G5A^ 
5D-\X,Y) . 5nX.Y)5G{V,W) _ 

6A^{z) '6G{V,W) 5A,{z) ^ > 

5A,{z) '5G{y,Wr'^^'^'^^ ' 

where in the last two lines the coordinates are shown explicitly (V and W are integrated). 
This is an integral equation that determines lis once the quantity 5S / 5G is known. Taking 
further functional derivatives of the SDE with respect to A^ (and setting A = results in 
a set of integral equations for n > 1 response functions. Note that S is directly obtained 
from r[G', A], so the SD and BS equations are fully consistent with each other. 

As we have seen, the combination of the SDE and the BSEs completely determine the full 
propagator and the response functions n(„_|_2)- The SDE is obtained from the 2PI CGEA, 
while the BSEs are obtained from the SDE by differentiation with respect to the external 
field. The important point is that, because the 2PI CGEA is invariant under external 
gauge transformations, both the full propagator and response functions obtained this way 
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are external gauge covariant. As we will show below, the EGI property of the 2PI CGEA 
implies that G and n(„+2)5 as obtained from T[G, A\, satisfy the WT hierarchy. This provides 
the required link between current conservation in nonlinear response and the external gauge 
invariance of the 2PI CGEA, and also a powerful and systematic way of studying nonlinear 
response in strongly interacting systems coupled to ideal reservoirs. 

To see the connection between EGI and the WT hierarchy, recall that external gauge 
invariance of the effective action means G[A'] = UG[A]U^^ . Inserting the expansion in 
powers of the external field given in Eq. (l38ll into both sides of this identity we get 



Z2 



'n P 

5^ ^ / . . . dnA^^{l) . . . A;jn)nf-,7(1, . . . , n; X, r) 

(56) 

UG[Q]u-^ + zA^f/n^f/-i + -A^A,f/nf f/-i + . . . 

I dl...dnA,,{l)...A,M[U{X)nl^-^^{l,...,n-X,Y)U-\Y)] . 

^ — n * 'J 



n=0 



Note that Eq. (156|) explicitly shows that the response functions n(„+2) are external gauge 
covariant. In particular, for an infinitesimal external gauge transformation U{X) ^ 1 + 
ieip{X) the transformed external field is 

A'^iX) = A,{X) + dMX) , (57) 

so Eq. (13^ becomes 



. ■ „ (58) 



n. . 

n=0 i=l 



where we have suppressed the arguments of II^^;^'^^" (1 . . .n\X,Y) and defined [dn] = dl . . .dn 
for brevity. Comparing terms of the same order in A^ on both sides of this expression we 
get 

U2iX,Y)+t [ rfi[9;^(i)]n^(x,r;i) = 

J (59) 

U2{X, Y) + ie['f{X)U2{X, Y) - U2{X, Y^Y)] 
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for the zeroth order term, and 

dW2[9^^(2)]nr(X,F;l,2)A^(l) = 

(60) 

J dle[v{X)Ui;{X, Y- 1) - n^(X, Y- 1)v{Y)]A,{1) 

for the hnear term. Higher order terms have a similar structure but can become quite in- 
volved. The main point to emphasize is that the EGI of the 2PI CGEA implies relationships 
among the response functions, shown explicitly in Eq. (1581) . 

Assuming that (p vanishes at infinity, we can integrate by parts the second term in Eq. 
f l59|) (zeroth order in A^) to obtain 

I j rfMi)9;n^(x,F;i) = 

- ^e[^(X)n2(X, Y) - n2(X, Y)^{Y)] 

which implies 

a^n^(x,y; z) = -en2(X,y)[5(X -z)- 5{Y - z)] (62) 

or, in more detail. 



(61) 



dl < T,3^{z)i,{X)i,{Y) >= 

e < T^iIj{X)^{Y) > [6{Y - z) - biX - z) 



(63) 



This is precisely the identity corresponding to n = 1 in the WT hierarchy given by Eq. ([T]), 
that is, Eq. ([7]). We see that, even at zeroth order in the external field, the relation between 
the three-point vertex and the two-point function, Eq. fl63l) . is satisfied due to the EGI of 
the 2PI EA. 

Similarly, for the linear term, Eq. fIBUl) implies the following identity 



rf292nr(X,F; 1,2)^(2) = 

(64) 

e[¥.(X)n^(X, Y- 1) - n^(X, Y- \)^{Y)\ 

which leads by the same calculation as in the zeroth order case to a WT identity between 
the three- and four-point response functions (equivalent to Eq. ([1]) for n = 2): 

diw:[x,Y-xz) = 

(65) 

e{b{J - z) - b{X - z)\W,{X,Y-\) . 



20 



It is clear that this procedure could be continued to higher order terms in A^, thus 
generating higher order WT identities. We note that the hierarchy obtained for the response 
functions n(„+2) is completely equivalent to that involving A(„), given in Eq. ([T]), as expected 
since, ultimately, they both enforce current conservation. As it can be seen from their 
definitions, Eqs. ([2]) and (1471) . the relation between vertex and response functions can be 
compactly written as 



nf-2^)"(X,F;zi...z„) = (-e)("-i)n l'"^ I)'^"«,^.)Af^)(X^2...^n,1^4---<;^i) (66) 

a=2 

for n > 2, and 

nf3)(X,F;^)=Af,)(X,F;^) for n = 1, ^^^^ 
n(2)(X, Y) = iG[0]{X, Y) for n = 0, 

where the operator D^°'{z'^, Za) is defined in Eq. ((31). Defining the operator (with n > 2 and 

2<a<n) 

n 

that appears in Eq. fl66l) acting on A|*^^, noting that this operator does not depend on zi 
(the "external" coordinate in A(„)), and using the property of the Green's function G„ 

({z;, zJ)G'„(l ...z^...n,l' ...z'^...n') = 
Fr, < T,^(l) . . . i^{z^) . . . ^(n)^(n') . . . ^P{z'^) . . . ^(1') > 

(69) 

=< T,fHz2) . ..f-{z^)i;{l)i;{l') >= 

nf^+2^)"(l5 ^2 • • • Zn) 

the equivalence between WT hierarchies can be proven using Eq. (166|) (the calculation is 
tedious but rather straightforward). 

Therefore, the external gauge invariance of the 2PI CGEA enforces the WT hierarchy 
necessary for current conservation beyond the expectation value level. By expanding the full 
propagator (solution to the SDE) in powers of the external field, we can calculate current- 
conserving response functions as solutions to the BSEs obtained from the SDE. The key 
point is that the EGI of the 2PI CGEA implies the covariance of the full propagator and 
the response vertex functions, and this results in the WT relationships among them. We 
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emphasize that such a procedure rehes on an expansion of G[A] in powers of the external 
field A, so it is valid only for weak external fields. 

We end up this section by indicating how to calculate the current induced by the external 
classical field. The average current which evolves by the action of the external field is given 
by Eq. (jS]), which can be re-expressed in terms of the lesser Green function G^^: 

< f{z) > = -e lim D^'iz, z')[A] < tlj\z')tlj{z) >= 

(70) 

ie lim D''{z,z')[A]G+4z,z') 

where /^'^[A] is given in Eq. ([5]) and one should recall that must be calculated in the 
presence of the external fielcl3. The lesser Green function appearing in the expression for 
the current can be calculated directly from the expansion given in Eq. ( |38l) . In this way, the 
induced average current can be systematically calculated to arbitrary order in the external 
classical field. The approach based on the 2PI EA and its loop truncation, adopted here, 
guarantees that the WT hierarchy is fulfilled. 



VI. APPROXIMATE 2PI EFFECTIVE ACTIONS AND CURRENT CONSERVA- 
TION 

As shown in Sec. Illlt the virtue of the 2PI CGEA method is that it provides a systematic 
way of encompassing interacting fields and current-conserving nonlinear response theory in 
a unified way. In particular, truncations to the loop expansion of the 2PI CGEA are, to 
arbitrary order, external gauge invariant. Therefore, any such approximation will preserve 
the WT hierarchy. In this section, we give a concrete example of an ad hoc approximation 
to the 2PI CGEA, i.e. not obtained from a truncation in the loop expansion, and analyze 
the consequences for current conservation. 

To be specific, we consider the following bare interaction vertex in Smt as appearing in 
Eq. ©: 

Uabcd = Cabcd\r^uU (71) 

where 

1 if a, 6, c, d = + 
Cahcd = S — 1 if a, 6, c, (i = — (72) 
otherwise 
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is a CTP tensor. In Eq. (17T!) we have grouped delta functions as follows 

S{rA- rB)S{rA- rc)S{rA- td) (73) 
= 6{Xa - XB)SiXA - Xc)6{Xa - Xn) 

and [a stands for spin projection) 

This vertex describes a "Hubbard-like" local interaction n|(r)nj^(r), where n^{r) is the elec- 
tron number operator at position r. 

It will be convenient for what follows to write the quantum correction to the 2PI EA as 

oo 

r2[G] = 5^r« (75) 

i=l 

where r2'' is of order f/'. We note that, due to the fact that fermion mean fields vanish, 
an expansion of r2 in powers of U is equivalent to a loop expansion (this is no longer true 
when symmetry breaking occurs, see e.g. Ref. Q). The 0{U) diagram corresponds to the 
double-bubble, which is two-loop, while the 0{U^) one corresponds to the basketball, which 
is three-loop. Both diagrams are given in Eq. ( l23l) . Since the CTP self-energy is obtained 
by functionally differentiating r2 with respect to the full propagator, the contributions to 
the self-energy can also be classified according to Eq. (1751) . According to Eq. fl7T]) . for i = 1 
we have 



^^AB -^^abi^A,XB) - -UabCdGcD 



=\ucabcd J2 J dXcdXn6{XA - Xb)6{Xa - Xc)6{Xa - Xn) x 

o"c.o"D (76) 

SaA<TB^<^Cf^DGcd{Xc, Xn) 
= -UCabcdS{XA — XB)Sa^aBGcd{XA, Xb) 

We note that the exchange part of T^ab is absent because the interaction occurs between 
particles with opposite spin projection. 
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From this equation, the (+, — ) components of the CTP self-energy are immediately ob- 
tained: 

SV'KXa, Xb) = U5{Xa - Xb)5 (TA<^B^++{-^Ai Xb) 
T.^11{Xa,Xb) = -U6{XA~XB)6.,,a^G-4XA,XB) (77) 
si^l = S« = 

where the last line follows from the definition of the CTP tensor Cabcd, and shows the well- 
known fact that the quasiparticle life-time is infinite in the Hartree approximation. 

One can proceed similarly for i = 2 to obtain (we omit the calculation since it is com- 
pletely analogous to the i = 1 case) 

e'^^1{Xa,Xb) = -iUH{Xa - XB)5,^.AG++{XA,XB)f 
eL'1(X^,Xb) = -iUH{Xa - XB)5,^.AG-{XA.XB)f 

(78) 

S (X^, Xb) = tUH{XA - Xb)6,^„, (X^, Xb)]' 
SL'|(X^, Xb) = iUH{Xa - XB)6,,,a^ [G-+(Xa, Xb)]' 

showing that this approximation includes fluctuation damping since the lesser and greater 
components of the self-energy are nonzero. The full CTP self-energy to order U"^ is then 
given 

As we have shown in previous sections, the truncation of the loop {U) expansion of 
'^2[G\ does not violate the external gauge invariance of the full 2PI EA, therefore provid- 
ing approximate propagator and vertex functions that satisfy the WT hierarchy. We will 
now discuss, in the context of the 2PI EA description presented in this work, an ad hoc 
approximation known in non-equilibrium perturbation theory applied to transport through 
quantum dots (see Refs. [llis^ and references therein), and show why it does not preserve 
current conservation. 

In the language of the 2P1 EA formalism, the approximation involves two separate steps. 
In the first one, the quantum correction r2 is approximated by the two- loop {0{U)) contri- 
bution, T^2^ . The 0{U^) self-energy is calculated self-consistently either from^^ the Hartree 
G^^ (lesser correlation) or by requiringii^ that the occupation of the central region evalu- 
ated with a renormalized S*^^) (but not the Hartree one) equals the one calculated (in the 
next step) from the propagator dressed with S^^-*. In the second step, S^^-* is calculated from 
Eq. ( ITHI) but with the full propagator Gab replaced by that calculated from S^^^ in the 
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first step. In simpler terms, in this approximation the internal lines in E^^^ are not the full 
propagator, as they are in the 2PI EA approach, but are either the Hartree or a similarly 
calculated propagator. It is easy to recognize that such approach can not, as it is, produce 
current conserving results in general, since it breaks the variational procedure by which the 
full propagator is obtained from a single (truncated) functional. 

The difference between the approximations to the 2PI EA based on a loop truncation and 
the ad hoc one can be best appreciated by comparing their SD equations. For the truncation 
to three loops we have 

Gab = D~/s + « [G] + zSj], [G] (79) 

where we have made explicit that S*^^'^) are functionals of the full CTP propagator Gab, 
given by Eqs. fl77|l and fl78|) respectively. For the ad hoc approximation we have instead the 
following system of self-consistent equations: 

9ab = Dab + ^^ab[9] + K 

~9AB = DAB + ^^%[9]+^^^Ik9] (80) 
9+~{r,t;r,t) = g+^{r,t;r,t) 

where g is the 0{U'^) propagator in this approximation scheme, and i^' is a source added 
in the calculation of g to enforce the condition expressed in the third line. This condition 
corresponds to requiring that the level occupation of the central region calculated from g 
and from g be equal. Note that if = then the 0{U) propagator, g, becomes the usual 
Hartree propagator. The equation of motion for g can be cast in a more familiar ioiwr^ 
using the properties of CTP propagators described in Sec. |Vl The result is (1 = (ri,ti)) 

^+_(1,2) = j dUA ^,,(l,3){S<^,+_(3,4) + sV'iM(3,4)}^,(4,2) 

^-+(1,2) = j rf3rf4^,(l,3){S^,_+(3,4) + sL'|M(3,4)}^,(4,2) (81) 

^,(1,2) = (7,(1,2) + j d3(i4^,(l,3)S(2)[(7](3,4)^,(4,2) , 

where the retarded 0{U'^) self-energy is 

(1, 2) = 0(1, 2) [S?i(l, 2) - Sjl(l, 2)] (82) 

and is given in Eq. (fT7|) . Using that the reservoirs are in equilibrium and non-interacting, 
and that the coupling matrix does not mix CTP branches, i.e. Tab = diag(T, T), the 
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components of in energy-momentum space can be written as 

= \T\'B+4p) = J2'^m\T\Mp)S{po-py2) 

j 

= \T\'B_M (83) 

= -J2'^m\T\'[l - nip)]5{po-py2) 
j 

where n{p) is Fermi-Dirac function, the momentum index j denote single-particle energy 
levels, and we have assumed that T{p) = T, independent of p. We note that the term 
describing the evolution of initial correlations is absent from Eq. flSTl) . since we are assuming 
that the system is in a non-equilibrium stationary (or periodic) state mantained by the 
external classical field. 

The equation of motion for the full propagator G has the same structure as the one for 
g, with the following replacements 

sJIM, sL'IM, ^ S_+[G], S,[G] (84) 

We will now consider the expectation value of the electronic current in the ad hoc ap- 
proximation discussed so far. As emphasized in Sec. [Tll the WT hierarchy given in Eq. ([1]) 
is a direct consequence of the equations of motion satisfied by the exact propagators. This 
can be shown quite simply by noting that the {n + 1)— point vertex function A|^^^, defined 
in Eq. ([2]), can be rewritten in terms of the operator D^^{z, z'), defined in Eq. ([S]): 

Af„)(l...n, l'...n';z) = 

- ei("+^) lim D^'iz, (1 ...nz,l'... n'z') 

The divergence of A|^^^ will therefore be given by the divergence of the operator D^{z,z'), 
which follows immedeately from its definition. The result is 



A{z) + A,(z')]G'(n+i)(l ...z;l'...z')}= (86) 
26 
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d,Al^il...n, l'...n';z) = 
^"elimC(^,^')G'(„+i)(l...^; l'...z') 



This equation shows that the divergence of the expectation value of the four-current j^{z), 
that corresponds to n = 0, is directly given by the action of the differential operator ({z, z') 
acting on {z, z') = G{z, z'). Using Heisenberg's equation of motion for the field operators 
ip, ip'^ it is easy to show^^'^^ that the exact propagator satisfy two differential equations 



[idt^ Ai{z)]G{z, z') = 6{z, j dl U{z, l)G2{zl, /!+) 

[~idt^, + - Ai{z')]G{z, z') = 5{z, z')-i j dl f/(l, z')G2{zl-, z'l) 
that, when substracted, result in 



(87) 



\imC{z,z')G{l...z]l' ...z') = Q , (88) 

whereby (9^ < j^{z) >= in the exact theory. 

In an approximate theory, the quantity —iU G2 appearing in Eq. (1871) is replaced by the 
product SG, and so the conservation of current at the expectation value is not immedeately 
guaranteed. It is well-known^>i^>i^>2^ that $— derivable approximations are such that < 
j^{z) >= 0. As we have shown in previous sections, the approximations to the 2PI EA 
based on truncations of the loop expansion guarantee not only this condition but the entire 
WT hierarchy (order by order in the expansion of G with respect to the external driving 
field). 

The approximation defined by Eq. (I80p can not be derived from a loop-truncated effective 
action, and hence Eq. (IHHj) is not expected to hold. To end up this section, we will briefly 
show how large is the divergence of the current expectation value in the ad hoc approxima- 
tion. From Eq. (IHTl) we have that, for an approximation based on a truncation of the loop 
expansion of the 2PI EA, the following identity holds: 

lim ({z, z')G{z, z') = 

y (89) 

/ dl {G'(z,l)S(2)[G'](l,2)-S(2)[G'](2,l)G'(l,2)} = 

where we note that the Hartree part of the self-energy, can be included in the left 

hand side of Eq. (IHTI) as an extra single-particle term, and, being local, its difference vanishes 
when z' z. Therefore, in Eq. fl89l) we only need to consider S'^^^ Putting G = g + G 
and S'^^^[G'] = + S in Eq. (189!) . and neglecting terms proportional to G and H, we 
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immedeately obtain 

d, < f{z) >= 

r (90) 

e j dl {~g{z,l)^^'^[9]{l,z)~E^'\g]{z,mi,z)} 

where j denotes the current calculated in the ad hoc approximation. Since g, S*^^^ [g] ~ 0(f/^), 
we conclude that the violation of current conservation in the mean is 0{U^) or higher. 

To conclude, we note that the combination of Eqs. (15^ and constitute a useful way 
of checking that a given approximation is conserving in the mean, and may provide some 
insight in the search of conserving approximations not based on loop-truncations. As we 
have shown, it also provides a way of calculating an upper bound to the violation of mean 
current conservation. 



VII. CONCLUSIONS 



In this work, we have determined the basic requirements that an approximation to a 
non-equilibrium many-body problem in an open and driven fermionic system must satisfy in 
order to achieve current conservation beyond the expectation value level. One of the most 
important results of this work is the close relation found between nonlinear response theory 
and the Ward-Takahashi hierarchy, necessary for current conservation. This connection, 
although already known^'^'^iiSiiiiiiii^^ was clearly displayed by using the closed-time path two- 
particle irreducible coarse-grained effective action to describe the central electrons coupled 
to ideal reservoirs. 

We have shown that the gauge invariance of the 2PI effective action with respect to trans- 
formations of the external classical field driving the system automatically implies the WT 
hierarchy among the propagator and vertex functions calculated from the effective action. 
More importantly for practical calculations, for every approximation to the 2PI effective 
action that results from truncating its loop expansion, the closed-time path propagator (ob- 
tained from the Schwinger- Dyson equation) and the the vertex functions (obtained from the 
Bethe-Salpeter equations) are such that the WT hierarchy holds. 

Using a simple example, we have also discussed, in the context of the 2PI effective action 
formalism, the relation between ad hoc approximations (not obtained from a truncation of 
the loop expansion) and current conservation. Using a general expression for the divergence 
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of the mean current, we have shown that in the ad hoc approximation considered current 
conservation at the expectation value level is violated at 0{U^) or even at higher order. 

In summary, we hope to have shown that closed-time path 2PI effective action techniques 
are a powerful and systematic method to study the nonlinear response in strongly correlated 
open systems in weak external fields, in a current-conserving way. We also hope to have 
shown the necessity, within the EA approach, of using loop-truncated approximations to 
the 2PI EA when going beyond the linear response regime. Our results may be of use in 
the theoretical study of quantum transport through interacting electronic pumping devices, 
which are nowadays receiving much attention. 
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Of course, no conflict arises for non-interacting particles, since in that case the lifetime of field 



22 for a discussion of second order response theory in non- 



excitations is infinite. See Ref. 
interacting systems. 

35 Wg refer to approaches to solve the problem that still belong to the two-particle irreducible 
class. Of course, one could in principle solve the problem of local symmetries (WT identities) 
by treating the dynamics of higher order correlations. This would require n-particle irreducible 



effective actions, with n> 2 (see Refs. 14 and 



26l ). However, we think that the relation between 



nonlinear response and current conservation is clearly appreciated in the 2PI formalism, which 
is also technically feasible as compared to n > 2 effective actions. It remains to be investigated 
the possibility of relating higher order response functions (and the response of higher order 
correlations) to higher order effective actions in the n-particle hierarchy. 

The causality of the response functions is easily proved in the physical representation, and 
follows directly from the structure of the generating functional of OTP correlation functions. 
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